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Abstract 

A mass and spring model is used to calculate the phonon mode dispersion for single wall carbon 
nanotubes (SWNTs) of arbitrary chirality. The calculated dispersions are used to determine the 
chirality dependence of the radial breathing phonon mode (RBM) density. Van Hove singularities, 
usually discussed in the context of the single particle electronic excitation spectrum, are found 
in the RBM density of states with distinct qualitative differences for zig zag, armchair and chiral 
SWNTs. The influence the phonon mode density has on the two phonon resonant Raman scattering 
cross-section is discussed. 

PACS numbers: 81.07.De,63.22.+m 
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Understanding how spatial confinement influences both a material's excitation spectrum 
and the density of excitation modes is imperative to accurately modeling the observable 
properties of the system. In particular, there has been much recent interest in the nanoscale 
confinement of phononsi. In one-dimensional (ID) single wall carbon nanotubes (SWNTs), 
there are various methods employed to calculate the phonon mode dispersion. Early work, 
following the successful calculations of the single electron bandstructure, projected the two- 
dimensional (2D) graphene phonon disperion relation to ID. It was quickly realized, though, 
that zone-folding did not capture all the phonon modes of SWNTs^. Following some ab- 
initio work^, a microscopic mass and spring model has been formulatecl^ to accurately 
capture the lattice dynamical properties of SWNTs. The mass and spring model permits 
calculation of both the phonon dispersion relation and phonon mode density for SWNTs of 
arbitrary chirality. 

In the optical characterization of SWNTs, one phonon resonant Raman scattering 
(IphRRS) is employed to measure a tube's electronic resonances, investigate various Ra- 
man active phonon modes and even determine the physical diameter and chirality of the 
tube under study^. Though IphRRS is a powerful tool for SWNTs characterization, con- 
servation of momentum allows only a single phonon to contribute to the measured signal 
at a particular scattered photon frequency. In contrast, in two phonon Raman scattering, 
multiple phonon pairs originating from the full phonon dispersion curve can contribute to 
the scattered signal. Therefore, the two phonon resonant Raman scattering (2phRRS) cross- 
section, unlike the IphRRS case, is proportional to the phonon joint density of states (DOS). 
The 2phRRS cross-section is nonzero for all scattered frequencies for which the product of 
the phonon joint DOS and transition matrix element is nonzero. Although 2phRRS has 
been discussed in the context of SWNTs , the discussion has only focused on calculating 
and little attention has been paid to the density of phonon modes^. It is clear then, 
in 2phRRS, an accurate model of both the phonon dispersion and joint DOS is important. 

To calculate the phonon dispersion and DOS we model the nanotube lattice as a collection 
of equal mass points, capable of displacing in all three spatial dimensions. We use Newton's 
law F = — /cx = ma to understand the dynamical evolution of each lattice site^. The 
governing equation is 

n'^Q = ]N{a,q;n,m)Q (1) 
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where (n, m) index the nanotube, a = 0,l,...,A^ — 1, A^= ^'^^n+m Xm+n) equals the number 
of graphene unit cells in the nanotube unit cell, q is the ID phonon crystal momentum 
along the tube axis, VL^ = muj'^/Ki is the rescaled phonon mode frequency where Ki is 
the force constant that characterizes the strength of the various couplings to be described 
below, Q = [Qa,p QA,e Qa,z Qb,p Qb,9 Qb,z] is the vector of fourier amplitudes for lattice site 
displacements with Qi j indicating the carbon atom at basis site i displaces in the direction 
j and the matrix dictating the dynamics is 

W = ^SiMutNN + S2M2ndNN + S-^Mrbb^ ■ (2) 

In W, MistNN is a 6x6 matrix characterizing the influence first nearest neighbor carbons 
atoms have in displacing a specific basis atom, M2ndNN is a 6x6 matrix characterizing the 
influence second nearest neighbor carbon atoms have in displacing a specific basis atom 
and Mrbb is a 6x6 matrix characterizing the restoring force on a given basis site resulting 
from bending the bond between neighboring carbon atoms. The above dynamics allow each 
carbon atom to displace in the p, 9, and z directions. In Eq. 0, Si determines the relative 
strength of the given coupling with respect to the coupling characterized by the spring 
constant Ki. Specifically, Si = ^ where Ki is the spring constant characterizing the i*'* 
restoring force. The parameter a serves as a band index for the phonons, where the number 
of phonon bands for a given nanotube is equal to three times the number of carbon atoms 
in the nanotube unit cell. For fixed a, the phonon mode dispersion relation is calculated by 
diagonalizing Eq. Q as a function of g as g ranges throughout the nanotube's first Brillouin 
zone. 

In addition to determining the dispersion relation, it is also possible to calculate the 
ID phonon density of states. Generally, the density of phonon modes is an important 
quantity when calculating transition rates for scattering processes mediated by phonons. We 
are particularly interested in the two phonon Raman scattering process. The two phonon 
mediated Raman scattering cross-section, assuming free electrons and holes, is^ 

7 2-ph 

= Cp%-of\2n,) . el; = - 2Q„ e,) \' (3) 

where we have assumed Stokes scattering, C is a constant, p^'^^g\2Vtp) is the two phonon 
density of states and [W^^j'^P is the transition probability from initial system state i, with 



3 



a single pump photon, to a final state /, with a single scattered photon and a pair of 
degenerate single phonons with opposite crystal momentum. We note, in Eq. Q our focus 
is on final states that contain two degenerate phonons of frequency Qp. In 2phRRS, it is also 
possible to have scattered photons at lOs = uji — {^pi ±^lp2), the sum or difference frequency 
of two nondegenerate phonons. In a 1-dimensional system, our expectation is the reduction 
of phonon phase space will hamper the likelihood of simultaneously satisfying both energy 
and momentum conservation in the sum or difference frequency generation process so we 
have ignored this possibility in Eq. (jS)). 

The utility in focusing on the degenerate phonon generation, is that we can approximate 
P%^s\'^^p) Pdc)s^\^p)^- With the single particle density of states defined as^ 

PDOs{^)dQ oc -^4^ (4) 

I dq I 

we can evaluate p^oos^i^p) f*-*^ ^ given phonon band with knowledge of a). In degen- 
erate 2phRRS, P^DOs^i^p) serves as an indicator of how many phonon pairs contribute to 
a specific scattered photon frequency. Specifically, the regions of zero slope in the phonon 
dispersion result in Van Hove singularities^ that enhance the 2phRRS cross-section. 

To illustrate the chirality dependence of the REM single phonon density of states, we 
study nanotube families 2n + m = 22 and 2n + m = 23 as well as subsets of both armchair 
and zig zag tubes. Our interest in family 22 stems from recent experimental work in 2phRRS 
from single SWNTs^^. In addition, family 22 and 23 yields a set of semiconducting SWNTs 
with chiral angles ranging from 0° to 30°. 

Before studying the radial breathing mode (RBM) DOS, we illustrate a typical phonon 
dispersion relation calculated from Eq. Figure ^a) is the full phonon dispersion relation 
for an (11,0) nanotube (family 22). The number of carbon atoms in the nanotube unit cell 
is 44 which results in 3 * 44 = 132 phonon modes. Also, the length of the primitive lattice 
vector along the tube axis is 3 times the carbon-carbon distance on the graphene lattice. To 
fix the frequency axis, we have followed Mahan^ and set Vt = ^/S equal to to = 1600 cm~^ 
(this fixes the value of the free parameter Ki). We also assumed in Eq. Q that si = 1, 
S2 = 0.06 and S3 = 0.024^. The calculated zone center frequency for the RBM of the (11,0) 
tube is 266 cm~^. The calculated value is in near perfect agreement with recently reported 
RBM-diameter relation of ojrbm = ^ + 18 where D is the nanotube diametei^, which 
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yields a zone center RBM frequency of 267.4 cm~^ for the (11,0) nanotube. 

Now we investigate the RBM DOS chirahty dependence. In Fig. ^b), we plot the 
RBM dispersion curve and RBM DOS for SWNTs of family 22. With the exception of 
the tube (11,0), which will be discussed when we focus on the zig zag tubes, the other 
tubes comprising family 22 have RBMs with mode densities concentrated entirely at the 
zone center. In addition, the chiral tubes of family 22 have RBMs that exhibit very low 
dispersion across the nanotube Brillouin zone. The Van Hove singularity in the RBM mode 
density of family 22's chiral tubes results in a 2phRRS cross-section that is entirely mediated 
by the zone center RBM. In this case, even though the momentum conservation has been 
relaxed in the two phonon scattering process, the phonon mode density prohibits phonons 
of appreciable crystal momentum from contributing to the 2phRRS cross-section. Family 
23 exhibits qualitatively identical behaviour to family 22. 

Figure 121 contains plots of the RBM dispersion, and its associated DOS, for semiconducting 
zig zag tubes (8,0), (10,0), (11,0) and (13,0) where we note that the RBM bandwidth 
decreases with decreasing tube diameter. There are two striking features in Fig. Efa) when 
compared to the chiral tubes of family 22. First, the bandwidth of the RBM is considerably 
larger than the chiral tubes, nearly 2 orders of magnitude in some cases. Second, zig zag 
tubes exhibit two sharp Van Hove singularities in their RBM mode density. If we define the 
zone boundary of the Brillouin zone as Qmax, we find that only wavevectors out to 0.03 * Qmax 
contribute to the zone center singularity. In 2phRRS, when the incoming laser frequency is 
commensurate with an electronic resonance, the two singularities in the RBM mode density 
should result in two peaks in the scattering cross-section, with a higher intensity peak 
associated with a frequency that is twice the zone center RBM frequency. 

Finally, we examine the RBM dispersion and DOS, plotted in Fig. Efb), for the (metallic) 
armchair tubes (7,7), (8,8), (9,9) and (10,10). In comparison to the previous zig zag tubes, 
the strongest Van Hove singularity for the RBM has shifted from zone center to the vicinity 
of the zone boundary. Approximately 0.01 * Qmax wavevectors contribute to the zone bound- 
ary singularity. The previous observation has an important consequence in 2phRRS. Not 
only will we observe a double peak in the scattering cross-section when the incoming laser 
frequency is tuned to an electronic resonance, but the more intense peak will be associated 
with a frequency that is equal to twice the zone boundary RBM phonon frequency. 

In summary, we have used a microscopic model of the lattice dynamical properties of 
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SWNTs to determine the chirality dependence of the RBM mode density. The importance 
of the phonon mode density for 2phRRS was discussed and illustrated. For chiral tubes of 
family 22 and 23 (with the exception of the (11,0) tube), we found narrow band RBMs with 
Van Hove singularities restricted to the zone center. In 2phRRS, we expect to see spectral 
features at scattered frequencies that are exactly 2 times the zone center RBM phonon 
frequency as determined by IphRRS. 

In contrast to the chiral tubes, the set of achiral SWNTs studied exhibited large band- 
width RBMs with 2 Van Hove singularities in their DOS. The zone center Van Hove sin- 
gularity dominated for zig zag tubes, whereas the zone boundary singularity dominated for 
armchair SWNTs. The appearance of two Van Hove singularities in achiral tube's mode 
density allows for the possibility of a 2 peak structure in the 2phRRS cross-section at fixed 
laser frequency. 
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FIG. 1: (a) The phonon dispersion relation for an (11,0) nanotube. The bold line is the RBM 
dispersion curve, (b) The RBM dispersion and RBM DOS for family 22 SWNTs. The number in 
the dispersion figures is the zone-center RBM frequency. In all figures, T is the magnitude of the 
primitive vector along the tube axis. 
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FIG. 2: (a) The RBM dispersion and RBM DOS for zig zag SWNTs (8,0), (10,0), (11,0) and 
(13,0). (b) The RBM dispersion and RBM DOS for armchair SWNTs (7,7), (8,8), (9,9) and 
(10,10). The number in the dispersion figures is the zone-center RBM frequency. In ah figures, T 
is the magnitude of the primitive vector along the tube axis. 
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